A 149  810 


AFWAL-TR-S4-4104 


\ 


PRELIMINARY  REVIEW  OF  STATISTICAL  TREATMENT  OF 
S-GLASS/ EPOXY  STRESS  RUPTURE  DATA 


MATERIALS  ENGINEERING 
UNIVERSITY  OF  DAYTON 
300  COLLEGE  PARK 
DAYTON,  OH  45469 


SEPTEMBER  1984 


»  • 


FINAL  REPORT  FOR  PERIOD  JUNE  1983  -  JUNE  1984 


i 


A 


1 

P 

< 


OJ 

o 

o 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED 


U  j 


C..« 

pn 


MATERIALS  LABORATORY 

AIR  FORCE  WRIGHT  AERONAUTICAL  LABORATORIES 

AIR  FORCE  SYSTEMS  COMMAND 

WRIGHT- PATTERSON  AIR  FORCE  BASE,  OHIO  45433 


DTIC 


NOTICE 


When  Government  drawings,  specifications ,  or  other  data  are  used  for  any  purpose 
other  than  in  connection  with  a  definitely  related  Government  procurement  operation, 
the  United  States  Government  thereby  incurs  no  responsibi lity  nor  any  obligation 
whatsoever ;  and  the  fact  that  the  government  may  have  formulated ,  furnished ,  or  in 
any  way  supplied  the  said  drawings ,  specifications ,  or  other  data,  is  not  to  be  re¬ 
garded  by  implication  or  otherwise  as  in  any  n.anner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights  or  permission  to  manufacture 
use,  or  sell  any  patented  invention  that  may  in  any  way  be  related  thereto. 


This  report  has  been  reviewed  by  the  Office  of  Public  Affairs  (ASD/PA)  and  is 
releasable  to  the  National  Technical  Information  Service  (NTIS) .  At  NTIS,  it  will 
be  available  to  the  general  public,  including  foreign  nations. 


This  technical  report  has  been  reviewed  and  is  approved  for  publication ■ 


,  7 

Aju. 


THEODORE  J.  REWHART,  Chief 
Materials  Engineering  Branch 
Systems  Support  Division 


BENNIE  COHEN,  Technical  Area  Manager- 
Materials  Integrity  Branch 
Systems  Support  Division 


FOR  THE  COMMANDER 


WARREN 
System 
Materi a 
,  .i  r  Force 


N ,  Chief 
port  Division 
Laboratory 
Wright  Aeronautical 


Laboratories 


"If  your  address  has  changed,  if  you  wish  tc  be  removed  from  our  mailing  list,  or 
if  the  addressee  is  no  longer  employed  by  you r  organization  please  notify  AFWA.l,/Hi,SF.  > 
W-PAL'3,  OH  45433  to  help  us  maintain  a  current  mailing  list". 

Copies  of  this  report  should  not  be  returned  unless  return  is  required  by  security 
considerations ,  contractual  obligations ,  or  notice  on  a  specific  document ■ 


UNCLASSIFIED 

StaiHiTv  I'lASl)!'  KlAllON  oi  TIMS  PAGE 


jl  ■  ijHlt  Y  C*  AS'olf  ICA1  I  ON 


REPORT  DOCUMENTATION  PAGE 


lb.  RESTRICTIVE  MARKINGS 


II  T  N  1.  I  ASS  1 1  ICA  I  l(‘N  Aul  HOUR  Y 


2l»  DCCl  AV»:»n  SA1ION  UOWNt.HADiNij  hCHl.OULL 


J  r\  HI  OHMIM,  ORGANIZATION  MhPOHl  NuMULPlS) 

AFWAL-TR-84-4104 


1.  nisi  Rl  E)UT  ION,  A  VAI  LABlU  I  Y  OI  REPORT 

Approved  for  Public  Release,  distribution 
uni  united. 


jb  MONITORING  OMGANI/A1  ion  REPORT  NUMBLR(S) 

APWAL-'I  U-H4-/.  I  Oh 


uj  NAMi  0‘  i'l  m  OUMiNij  Orti.AN) NATION 

Materials  Engineering 
School  of  Engineering 


t't  AODHfc  SS  (i  it  \ .  .S  fiifr  and  ///*  ('ode  l 

University  of  Dayton 
300  College  Park 
Dayton,  Oil  46469 


fa  NAM tOf  MONITORING  ORGANIZATION 


APWAI7M1.se 


7h  ADDRESS  H'lty.  S tat*'  and  ///'  Cud ef 

WPAPH ,  Oi!  4V<33-(o3J 


8b  OFF  ICE  S  MBOL  9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
ilf  applicable)  r^^rir  r-  r  o 's  r» 


tk  AOL'HI  SS  .(  ./•  .  slab  and  *.IV  C>de 


E  3  3  6 1 5  -  B2  -  C  -  50  3  9 


10  SOURCE  Of-  FUNDING  NOS 


(AISC) 

Wri ght-Patterson  A f B ,  0;l  45433 


11  1  i  T  i  t  Itu  mitt-  Vcu'ifv  (’luixi/ituliiiM 

PRELIMINARY  REVIEW  OF  STATIST  I  CAL  TREATMENT  OF 


12  I’L  ITS *JN  A l  AOtiinmS) 

Gerald  Shaughnessy  and  James  A.  Snide _ _ _ 

1 3.1  nLI’OUT  |  13b  tTmV  cov:.  a-'.)  I  14  DA1I-  nf  REPORT  Itr.  ,ll„, 

Final  Report  |  J_?8 3, 0 June  84 1  September  1984 


PROGRAM 

PROJECT 

TASK 

ELEMENT  NO 

NO. 

NO. 

93046 

15 

WORK  UNIT 
NO 


14.  DAI  I-  Of-  REPORT  »  >  r  .  ,\u>  ,  Day  /  1  b.  I'AUt  LUUiMl 

September  1984  13 


16  supplementary  notation 

Please  see  other  side  of  this  page  for  continuation  of  the  title. 


COSAT  i  COOLS 


Id  SUBJECT  TERMS  i  (\j/i  fume  on  if  necessary  and  identify  by  block  number  I 

''Composites.’  S-GlassL 

Stress  Rupturec  Epoxy, 

Statistical  Analysis',  ^ 


VJ  a'us  t  ^Ai'l  -t  i‘’>nm  <*n  fi’iTrs.-  if  »n’.  **.saiirv  unit  menfiiy  by  hiiwlt  number  > 

'>The  purpose  of  this  report  i:.  to  provide  statistically  derived  tabular  design  data  to 
pi  edict  tne  estimated  stress  rupture  life  of  unidirectional  E-glass/epoxy  composites 
subjected  to  sustained  tensile  loading.  The  maximum  likelihood  estimated  life  (MLEL)  is 
tabulated  for  various  stress  levels  of  the  fiber  and  quantile  level.  The  lower  bound 
lifetime  with  confidence  intervals  of  95  and  99  percent  are  also  tabulated  for  the  various 
tibei  sti ess  and  quantile  levels.  The  estimated  lifetimes  were  calculated  from 
accelerated  stress  rupture  tests  uf  7.  T,  Chiao  (Lawrence  Livermore  Laboratory)  who  used 
uniaxial  strand  samples;  therefore,  the  estimated  lifetimes  reported  are  considered  to  be 
conservative.  1  .  L\x,  ^  l_ 


20 

UN 

DIS1  UUiUI  K>N  A  V  All  AUl  l  '  1  t  OI  A13SIRACT  j 

CL  ASSII  'L  D  UNl  iMl  T  [  D  I  J  SAMt  A$  MPT  U  DllCUSbRS  O  ! 

_ _ _ _ _ j 

21  Al’ST  RACT  SL  GUI  1 1  S  CLASSlFU 

CATION  1 

17  a 

N  A  Ml  OI  Ml  S  1*0  NS  lull  INfU  VI  G'.JAl. 

221'  TELEPHONE  NUMBER 
l/nr/ui/*  1  ’fu  ■ 

(513)  239-3241 

22c  Off  ICE  SYMBOL 

KL  262 

DD  rORM  1473,  83  APR 


EDITION  OF  1  jAN  73  IS  OBSOl  ETE 

(OVER) 


UNCLASSIFIED  _ 

$E  CUR  l  T  Y  CL  ASS  IF  iCATljN  OF  T  HlS  PAGE 


UNCLASS  1 1  1  LI) 

SECURITY  CLASSIFICATION  of  this  page  ^ ___ ___  ___ 

Item  No.  11.  Title  -  S-GLASS/EPOXY  STRC3G  RUPTURE  C.‘lTA 


PKEFACt 


lliis  report  coverb  work  performed  under  Air  Force  Contract  No. 
F33blb-8P-C~S03y.  The  work  was  administered  under  the  direction  of  the 
Materials  Laboratories ,  Air  Force  Systems  Command,  Air  Force  Wright 
Aeronautical  Laboratori es ,  Wri ght-Patterson  Air  Force  Base,  Uhio. 

Mr.  Bennie  Cohen  (AFWAL/MLSA)  was  the  Program  Project,  engineer. 

The  Principal  Investigator  on  this  program  was  Ur.  James  A.  Snide, 
Director  of  Graduate  Materials  Engineering  at  the  University  of  Dayton. 


DTIC 


ELECTE 

JAN  2  2  1985 


B 


Accession  For 

NTIS  GRA&I 
DTIC  T :.i  □ 

Unfunicva'.ced  0 

Jus  1 1  i  i  e  .■  1 1  i  on_. — - - - 


By - 

D 1  a  t  r  1 1. 1 1 1  o  i  _ _ 

Availability  Codes 
Avail  and/or 
Dlst  1  Spaciul 


/IT 


TABIC  or  CONTENTS 


I.  INTRODUCTION . 1 

II.  BACKGROUND . ? 

III.  ANALYSIS .  4 

IV.  RESULTS .  8 

V.  CONCLUSIONS .  12 

REFERENCES .  13 


v 


LIST  OF  ILLUSTRATIONS 


I  1GUUL  PAGE 

1.  Graphs  of  the  Wei  bull  Probability  Density  Functions  ( L  -  1 )  ....  5 

LIST  OF  TABLES 

TABLE  PAUL 

1.  SUMMARY  OF  STRESS  RUPTURE  DATA  FOR  S-GLASS/EPOXY  STRANDS .  3 

2.  ESI [MATED  LIFETIME  OF  S-GLASS/EPOXY  COMPOSITE  UNDER  SUSTAINED . 

TENSILE  LOADING . y 

3.  ESTIMATED  LIFETIME  OF  S-GLASS/EPOXY  COMPOSITE  UNDER  SUSTAINED . 

TENSILE  LOADING.... . 1U 

4  ESTIMATED  LIFETIME  OF  S-GLASS/EPOXY  COMPOSITE  UNDER  SUSTAINED . 

TENSILE  LOADING . 11 


INTRODUCTION 


In  May  1981,  an  aircraft  sustained  major  damage  when  a  high-pressure 
nitrogen  storage  bottle  installed  in  the  aircraft  ruptured  without  warning. 
The  bottle  was  a  filament-wound  S-gl ass/epoxy  composite  structure  containing 
nitrogen  at  approximately  5000  psi  pressure.  The  pressurized  bottle  is  used 
to  inflate  the  aircraft  crew  module  impact  attenuation  bag  and  the  pressure  in 
the  bottle  is  not  released  during  the  service  life  of  the  bottle  unless  the 
attenuation  bag  is  employed.  The  bottle  that  tailed  had  been  in  service  and, 
therefore,  pressurized  for  eight  years  when  tne  failure  occurred.  The 
Accident  Investigation  board  at  Sacramento  Air  Logistic  Center  concluded  that 
the  failure  was  due  to  a  stress  rupture  failure  of  the  nitrogen  bottle. 

The  stress  rupture  of  composites  under  sustained  tensile  loading  has  been 
extensively  studied  by  T.  T.  Chiao  and  associates  (References  1-2)  at  the 
Lawrence  Livermore  National  Laboratories  during  the  past  decade.  This  work 
has  shown  that  the  stress  rupture  life  decreased  as  the  stress  developed  in 
the  fibers  is  increased  relative  to  the  ultimate  strength  of  the  fibers.  For 
example,  the  stress  rupture  test  of  one  hundred  S-glass/epoxy  composites 
loaded  to  83.5%  of  the  ultimate  strength  of  the  5-glass  fiber  failed  in  the 
time  range  of  0.01  hours  to  8=55  hours.  A  similar  group  of  Strand  samples 
loaded  to  65%  of  the  fiber  strength  exhibited  initial  failure  at  8.6  nours  and 
all  the  samples  had  failed  at  5014  hours.  The  fibers  in  the  tailed  bottles 
were  loaded  to  36-40%  of  their  ultimate  flow  strength;  therefore,  there  was  a 
3-4%  probability  of  stress  rupture  failure  of  the  bottle  after  eight  years  of 
service.  As  a  result  of  this  investigation,  all  the  remaining  bottles  of 
similar  service  were  replaced  and  the  old  bottles  were  destroyed. 

The  Sacramento  Air  Logistic  Center  recommended  that  the  appropriate  AKSC 
design  handbooks  should  be  updated  to  provide  information  concerning  the 
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stress  rupture  failure  composites  subject  to  sustained  tensile  loading. 
This  report  is  an  initial  attempt  to  summarize  the  experimental  stress  rupture 
data  on  S-gl ass/epoxy  of  T.  T.  Chiao  into  tabular  form  which  can  then  be 
incorporated  into  the  appropriate  AFSC  design  manuals.  In  addition,  the 
methodology  of  the  statistical  analysis  of  strand  data  to  calculate  the 
Maximum  Likelihood  Estimate  Life  (MLEL)  and  lower  bounds  for  various  fiber 
stress  levels  is  discussed. 

11.  BACKGROUND 

In  the  report  (Reference  I)  by  T.  T.  thiao,  R.  E.  Glaser  and  K.  L.  Moore, 
approxi mately  10  years  of  stress-rupture  data  for  impregnated  fiber  strands 
were  stati stical ly  analyzed.  Tne  data  were  collected  at  seven  stress  levels 
ranging  from  480  ksi  to  189  ksi  representing  values  from  84%  to  33%  of  strand 
failure  stress.  The  S-gl ass/enoxy  strand  data  that,  are  checked  in  the  above 
report  are  listed  in  Table  1. 

Any  statistical  treatment  of  the  experimental  data  must  consider  the 
fol lowing  di fficulties: 

1.  The  data  were  collected  at  seven  different  stress  levels  (o)  and  the 
parameters  (a  and  p)  of  the  statistical  model  for  the  strand  lifetime 
are  dependent  on  the  stress. 

2.  The  nature  of  the  relationship  between  o  and  the 
parameters  {a  aid  p)  of  the  strand  lifetime  distribution  must  be 
determined  from  the  data. 

3.  There  is  little  data  for  values  of  o  below  b0%  UT$. 

The  probability  distribution  of  the  time  to  failure  for  a  strand  at  a 
given  stress  level  is  determined  if  u  and  p  are  known.  The  usual  method  of 
analysis  is  to  estimate  a  and  p  separately  for  each  stress  level.  The 
difficulty  of  using  this  approach  for  the  S-glass  data  is  that  very  little 
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TABLE  1 


summary  of  stress  rupture  data 

FUR  S-GLASS/EPOXY  STRANDS 


SIKLSS,  KSI 

PERCENTILE 

UF  FAILURE 
STRESS  (%  UTS) 

#  STRANDS 
TESTED 

#  EXACT 
FAILURE 
TIMES 

#  GROUPED 
IA1LURE 
TIMES 

^CENSORED 

times 

479.5 

(84) 

117 

114 

3 

0 

426.3 

(75) 

104 

104 

0 

U 

373.0 

(65) 

115 

105 

0 

10 

332.0 

(58) 

98 

21 

13 

64 

286.0 

(50) 

46 

4 

3 

39 

228.8 

(40) 

100 

2 

0 

98 

188  =  7 

( 

v  •*/ 

100 

1 

J- 

A 

V 

99 

data  arc  available  tor  stress  levels  below  28b  k s i ,  and  no  usetul  statistical 
estimates  can  be  constructed  using  only  the  data  tor  the  stress  levels  ot 
interest  below  286  ksi. 

What  is  new  in  the  statistical  treatment  of  this  data  is  that  a  model  is 
used  that  incorporates  the  dependence  of  the  parameters  a  and  ft  on  the  stress 
level  in  the  Weibull  time  to  failure  distribution.  Instead  of  estimating  the 
parameters  «  and  ft  separately  at  each  stress  level,  new  parameters  are 
introduced  that  model  the  dependence  of  e;  and  ft  on  the  applied  stress  level. 
Those  new  parameters  are  then  estimated  from  all  the  data  using  the  method  of 
maximum  likelihood.  Once  the  parameters  of  the  model  relating  the  dependence 
of  a  and  ft  on  the  stress  level  have  been  estimated  from  the  data,  estimates 
of  quantile  points  in  the  failure  time  distribution  can  be  determined. 

IIL  ANALYSIS 

The  resulting  lifetimes  of  the  stress-rupture  tests  at  each  given  stress 
level  were  modeled  as  a  Weibull  distribution.  The  distribution  function  for  a 
Wei bul 1  distribution  is 

F(t)  =  1  -  exp[-(t/ft)'x]  (1) 

where  f-(t)  is  the  probability  that  the  lifetime  of  a  given  strand  subject  to  a 
stress  (o)  will  be  less  than  or  equal  to  t,  ana  u  ana  ft  are  the  shape  and 
scale  parameters  of  the  Weibull  model  to  be  estimated  from  the  experimental 
data.  These  parameters,  a  and  ft  of  the  Weibull  model,  vary  with  the  given 
stress  level  (o). 

An  insight  into  the  meaning  of  the  parameters  a  and  ft  in  the  Weibull 

distribution  of  the  S-glass  data  can  be  obtained  by  noting  that  the  mean 

lifetime  under  a  Weibull  model  is  pr(l  +  -) ,  where  I'(l  +1)  can  be  found  from 

a  <jr 

tables  of  the  gamma  function  (Reference  3).  For  the  8-glass  data,  u  was 
approximately  one  for  all  stress  levels.  The  mean  life  of  a  strand  is 
therefore  approxi mately  ftT(2)  =  ft. 
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The  reliability  function  R(t),  which  is  the  probability  that  the  lifetime 
of  a  strand  subject  to  a  stress  (o)  will  exceed  t,  is  related  to  the 


distribution  function  as  follows: 

R(t)  =  1  -  F (t )  =  exp  [-  (i)^  U) 

Once  the  parameters  a  and  b  are  estimated,  the  reliability  can  be 
determined  for  any  stress  level  o. 

To  illustrate  the  effect  of  the  a  parameter  on  the  shape  of  the  Weibull 
lifetime  distribution.  Figure  1  below  shows  how  the  shape  of  the  Weibull 
density  changes  from  an  exponential  type  curve  when  a  is  1  to  an  almost  bell¬ 
shaped  curve  for  a  =  4. 


Jxu) 


Figure  1.  Graphs  of  the  Weibull  probability  density  functions  (u  =  1). 

The  basic  method  of  analysis  used  in  this  paper  to  estimate  the  values 
of  a  and  B  of  the  Weibull  model  is  known  as  the  method  of  maximum 
likelihood.  In  this  method,  the  joint  distribution  of  the  observed  data 
(called  the  likelihood  function)  is  regarded  3s  a  function  of  the 


parameters  a  and  p.  The  values  a  arid  B  that  maximized  the  likelihood  function 
are  referred  to  as  maximum  likelihood  estimates  of  the  parameters  a  and  p. 


Those  values  can  then  be  substituted  into  the  reliability  equations  to  create 
the  tables. 

As  stated  earlier,  the  parameters  a  and  u  vary  with  each  stress  level. 
The  following  two  relationships  were  used  to  describe  how  a  and  p  varied  with 


the  stress  level  for  the  data  being  analyzed. 


-r=~r  =  0.  +  O9o,  (3) 

a(o)  1  i  '  ' 

1 n ( 3 { o ) )  =  03  +  O4o  +  05a2  +  O^o3  (4) 

On  are  constant  coefficients  that  are  to  be  determined.  Other 
representations  of  the  relationship  between  a  and  3  and  the  stress 
level  (o)  were  considered  but  the  linear  model  for  a  and  cubic  model  for  3  fit 
best  with  the  experimental  data  (Reference  4). 

furthermore,  to  insure  that  8(0}  is  a  decreasing  function  of  0,  the 
constraint  0^  =  -300(0^  +  225o6)  was  imposed.  This  constraint  is  obtained  by 
taking  the  derivative  of  Eq.  4  and  setting  it  equal  to  zero  when  0  =  150, 

making  f?(o)  a  decreasing  function  of  a  for  values  of  o  from  500  ksi  to  150 

ksi.  This  constraint  is  derived  as  follows, 
b(u)  =  1 n[p(o )] 

2  3 

b(o)  =  +  O^o  +  O^o  +  G^o 

and  taking  the  derivative 

b‘(o)  =  0^  +  20^o  +  30^0 3 

and  setting  it  equal  to  zero  while  0  =  150  ,o4  is  obtained. 

04  -  -300(O5  +  22506)  (5) 

Using  these  representations  of  a(a)  and  8(0)  ,  Eqs.  3  and  4  were  then 

substituted  into  the  likelihood  equation  which  was  then  maximized  as  a 


function  of  the  six  variables  (Oj,  C,,,  C>3,  c>4,  0&,  og).  The  maximization  was 
subject  to  the  constraint  of  Eq.  5  (Reference  3). 
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The  maximum  likelihood  estimates  (o^,  02»  0^,  G^,  Og,  Og)  are  then 
obtained  by  setting  five  partial  derivatives  of  the  likelihood  equations  equal 
to  zero  and  solving  the  resulting  five  equations  and  five  unknowns.  (There 
are  only  five  unknowns  since  is  a  function  of  o 5  and  ©g.) 

This  system  of  equations  that  must  be  solved  is  nonlinear  and  presents  a 
difficult  numerical  analysis  problem.  R.  E.  Glaser,  developed  a  Fortran 
program  to  so1 ve  this  system  of  equations.  However,  his  method  is  highly 
dependent  upon  close  initial  estimates  of  (0^,  02,  Og,  0^,  Og,  Og).  If  good 
initial  estimates  are  not  available,  the  program  may  not  converge  to  the 
maximum  vector  quantities.  To  obtain  good  initial  estimates  of  the  six 
unknown  parameters  the  following  technique  is  used.  Estimates  of  a  and  p  are 
obtained  for  each  different  stress  level  by  the  maximum  likelihood  method. 
Then  first  and  third  degree  polynomials  are  fit  to  the  a  -  a  and  0  -  P  pairs, 

A  /V 

respectively,  (values  of  a  and  e  obtained  at  each  stress  level)  by  the  method 
of  least  squares.  The  least  square  estimations  of  the  coefficients  in  the 
linear  and  cubic  models  are  then  used  as  initial  estimates  of  the  unknown 
parameters  (0^,  0^,  Og,  0^,  Og,  ©g)  in  the  Fortran  program  developed  by 
R.  E.  Glaser. 

Once  the  parameters  a  and  8  of  the  Wei  bull  model  have  been  estimated  by 
the  above  procedure,  quantile  points  in  the  lifetime  distribution  of  the 
S-Gl ass/Epoxy  composite  can  be  estimated.  A  quantile  is  a  point  in  the 
distribution  that  has  a  given  percentage  of  the  distribution  less  than  it. 
For  example,  the  0.1  quantile  in  a  lifetime  distribution  is  a  time  such  that 
the  probability  of  a  strand  failing  before  that  time  is  0.1.  Usiny  the 
maximum  likelihood  estimates  («  and  8),  the  pth  quantile  in  the  lifetime 
distribution  is  obtained  by  setting 
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p  =  1  -  exp  [-  {^)  ]  (b) 

and  solving  for  t.  The  solution,  yielding  the  pth  quantile  is: 

t  =  exp  [~  1 n( -1 n( 1  -  p))  +  ln(B)]  (7) 

a 

IV.  RESULTS 

Table  2  was  generated  from  Eq.  7  and  gives  the  maximum  likelihood 
estimated  life  (MLEL)  for  failure  probabi 1 i ti es  ranging  from  0.0001  to  0.1  and 
for  stress  levels  ranging  from  26%  to  45%  of  the  ultimate  flow  stress  (UTS). 

For  example,  at  a  stress  level  of  30%  (UTS)  and  a  quantile  value  of  .001, 
the  table  value  is  20.49  years.  This  means  that  the  maximum  likelihood 
estimate  of  the  point  in  the  distribution  of  lifetimes  under  which  there  is  a 
probability  of  failure  of  0.001  is  20. 49  years.  This  estimate  has  uncertainty 
in  it  since  it  is  based  on  estimates  of  a  and  u  obtained  from  the  experimental 
data.  Tables  3  and  4  give  lower  confidence  bounds  at  95%  and  99%  levels  of 
confidence,  respectively,  for  quantiles  in  the  lifetime  distribution  under 
which  the  probability  of  a  strand  failing  is  either  0.0001,  0.001,  0.01,  or 
0.1.  The  second  set  of  tables  are  more  conservative  and  reflect  the  size  of 
possible  errors  in  the  estimation  procedures.  At  the  95%  level  of  confidence, 
the  values  in  the  table  were  constructed  so  that  there  is  a  probability  of 
0.95  that  the  true  quantile  points  in  the  lifetime  distribution  are  not  lower 
than  the  table  values.  To  illustrate  further  the  use  of  the  tables,  suppose 
that  an  estimate  of  the  0.001  quantile  point  in  the  lifetime  distribution  of  a 
strand  subject  to  a  stress  of  30%  (UTS)  is  desired.  This  is  a  point  in  the 
distribution  below  which  there  is  a  probability  of  0.001  that  a  strand  will 
fail.  Using  Table  2  for  the  MLEL,  the  estimate  is  found  to  be  20.49  years. 
This  is  the  maximum  likelihood  estimate  of  the  0.001  quantile  of  this 
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TABLE  2 


ESITMATED  LIFETIME  OF  S-GLASS/EPOXY  COMPOSITE 
UNDER  SUSTAINED  TENSILE  LOADING 

MLEL  CALCULATIONS  (YEARS) 

(MAXIMUM  LIKELIHOOD  ESTIMATED  LIFE) 


QUANTILE  LEVEL 


%  Fiber 
Stress 

Stress 
( KS I ) 

0.0001 

MLEL  ( YRS) 

0.001 

MLEL  (YRS) 

0.01 

MLEL  (YRS) 

0.1 

MLEL  (YRS) 

26 

148.7 

1.87 

26.49 

375.8 

5600.8 

27 

154.4 

1.82 

25.86 

368.9 

5528.6 

28 

160.2 

1.72 

24.57 

352.4 

5310,1 

29 

165.9 

1.58 

22.73 

327.8 

4966.6 

30 

171.6 

1.42 

20.49 

297.1 

4527.0 

31 

177.3 

1.24 

18.02 

262.7 

4024.5 

32 

183.0 

1.06 

15.46 

226.7 

3492.3 

33 

188.8 

0.88 

12.97 

191.1 

2960.4 

34 

194.5 

0.72 

10.63 

157.5 

2453,5 

35 

200.2 

0.57 

8.52 

127.0 

1989.5 

36 

205.9 

0.45 

6.69 

100.3 

1579.7 

37 

211.6 

0.34 

5.15 

77.6 

1229.2 

38 

217.4 

0.26 

3.89 

58.9 

938.1 

39 

223.1 

0.19 

2.88 

43.9 

702.7 

40 

228,8 

0.14 

2.10 

32.1 

517.1 

41 

234.5 

0.10 

1.50 

23.1 

374.1 

42 

240.2 

0.07 

1.06 

16.4 

266.2 

43 

246.0 

0.05 

0.73 

11.4 

186.6 

44 

251.7 

0.03 

0.50 

7.8 

128.9 

45 

257.4 

0.02 

0.34 

5.3 

87.8 
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TABLE  3 


ESTIMATED  LIFETIME  OF  S-GLASS/EPOXY  COMPOSITE 
UNDER  SUSTAINED  TENSILE  LOADING 

LOWER  BOUND  CALCULATIONS  (YEARS) 

95%  CONFIDENCE  INTERVAL 


QUANTILE  LEVEL 


%  FIBER 

STRESS 

0.0001 

0.001 

0.01 

0.1 

t 

STRESS 

(KSI) 

LB  (YRS) 

LB  (YRS) 

LB  (YRS) 

LB  (YRS) 

;  '.-'j 

26 

148.7 

0. 

27 

154.4 

0. 

28 

160.2 

0. 

29 

165.9 

0. 

30 

171.6 

0. 

31 

177.3 

0. 

32 

183.0 

0. 

33 

188.8 

0. 

34 

194.5 

0. 

35 

200.2 

0. 

36 

205.9 

0. 

37 

211.6 

0. 

36 

217.4 

0. 

39 

223.1 

40 

228.8 

41 

234.5 

42 

240.2 

0. 

43 

246.0 

44 

251.7 

0. 

45 

257.4 

1991.0 

1972.7 
1906.4 
1798.2 

1656.7 


1133.1 


TABLE  4 


ESTIMATED  LIFETIME  OF  S-GLA3S/EP0XY  COMPOSITE 
UNDER  SUSTAINED  TENSILE  LOADING 

LOWER  BOUND  CALCULATIONS  (YEARS) 

99%  CONFIDENCE  INTERVAL 


QUANTILE  LEVEL 


%  FIBER 
STRESS 

STRESS 

(KSI) 

0.0001 

LB  ( YRS) 

0.001 

LB  (YRS) 

0.01 

LB  (YRS) 

0.1 

LB  (YRS) 

26 

148.7 

0.05 

1.7 

53.7 

1297.5 

27 

154.4 

0.05 

1.8 

54.2 

1287.6 

28 

160.2 

0.05 

1.8 

53.3 

1247.5 

29 

165.9 

0.05 

1.7 

51.1 

1180.8 

30 

171.6 

0.05 

1.7 

47.8 

1092.7 

31 

177.3 

0.05 

1.5 

43.7 

989.3 

32 

183.0 

0.04 

1.4 

39.0 

876.9 

33 

188.8 

0.04 

1.2 

34.1 

761.4 

34 

194.5 

0.03 

1.1 

29.2 

648.1 

35 

200.2 

0.03 

0.9 

24.4 

541.1 

36 

205.9 

0.02 

0.7 

20.1 

443.4 

37 

211.6 

0.02 

0.6 

16.2 

356.9 

38 

217.4 

0.02 

0.5 

12.8 

282.3 

39 

223.1 

0.01 

0.4 

10.0 

219.5 

40 

228.8 

0.01 

0.3 

7.6 

168.0 

41 

234.5 

0.01 

0.2 

5.7 

126.5 

42 

240,2 

0.01 

0.2 

4.2 

93.9 

43 

246.0 

0.00* 

0.1 

3.1 

68.7 

44 

251.7 

0.00* 

0.1 

2.2 

49.5 

45 

257.4 

- rce  TUi 

0.00* 

ITKi  ~h d  tni i hi  Tn 

0.1 

1.6 

35.3 

1' 


distribution.  Due  to  sampling  error,  this  estimate  may  be  somewhat  higher  or 
lower  than  the  true  0.001  quantile.  Using  the  second  set  of  tables  (Table  4), 
it  is  found  that  with  99%  confidence  the  0.001  quantile  is  not  below  1*7 
years. 

V.  CONCLUSIONS 

In  using  the  tables  developed  in  this  review  for  design  of  pressure 
vessels  for  Air  Force  applications,  it  should  be  kept  in  mind  that  the 
statistical  analysis  that  resulted  in  Tables  2,  3  and  4,  are  based  solely  on 
strand  data.  Various  arguments  have  been  presented  relative  to  the  use  of 
strand  data  to  predict  the  behavior  of  filament-wound  pressure  vessels. 
Strand  data  are,  however,  the  only  extensive  data  available  on  stress  rupture 
of  S-Gl ass./Fcoxy  composites  and  should  be  a  consideration  in  the  design  of 
pressure  vessels. 
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